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A previously proposed variational approach for momentum-conserving systems [J. Liu et.al., Phys. Rev. 

E 91, 042910 (2015)] is extended to systematically investigate general momentum-nonconserving nonlinear 
lattices. Two intrinsic identities characterizing optimal reference systems are revealed, which enables us to 
derive explicit expressions for optimal variational parameters. The resulting optimal harmonic reference systems 
provide information for the band gap as well as the dispersion of renormalized phonons in nonlinear lattices. 
As a demonstration, we consider the one-dimensional 0 4 lattice. By combining the transfer integral operator 
method, we show that the phonon band gap endows a simple power-law temperature dependence in the weak 
stochasticity regime where predicted dispersion is reliable by comparing with numerical results. In addition, an 
exact relation between ensemble averages of the </> 4 lattice in the whole temperature range is found, regardless 
of the existence of the strong stochasticity threshold. 
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I. INTRODUCTION 

The understanding of the microscopic foundation of 
Fourier’s law is a longstanding problem in statistical physics. 
One of the promising ways along this direction is the study 
of heat conduction in low dimensional nonlinear lattices [1- 
3]. Among the investigations, an interesting perspective is to 
find a nonlinear lattice which displays normal heat conduc¬ 
tion, namely, where Fourier’s law holds. After numerous nu¬ 
merical simulations [4—9], a consensus that one-dimensional 
(ID) nonlinear lattices with on-site potentials can establish 
normal heat conduction has been reached. Unraveling the mi¬ 
croscopic mechanism of heat transfer across such momentum- 
nonconserving (MNC) lattices thus may provide crucial in¬ 
sights to the foundation of Fourier’s law. 

Phenomenologically [10], the heat conductivity k can be 
approximately obtained via the Peierls formula, i.e., k = 
'^2 k CkVklk , here k is the phonon wave number, the spe¬ 
cific heat of the phonon mode, Vk its phonon group veloc¬ 
ity and Ik the corresponding mean free path (MFP). Follow¬ 
ing this spirit, a theoretical proposal for phonon heat transport 
based on renormalized phonon has been carried out [11, 12]. 
Nevertheless, flaws in its phonon description [13] may pro¬ 
duce wrong phonon information, which limits its ability and 
generality. Trying to understand the process of heat conduc¬ 
tion from the phonon point of view, theoretical investigations 
on the phonon properties seems to be a necessary step. Thus 
a unified phonon theory which is capable of dealing with var¬ 
ious nonlinear lattices is of primary interest and importance 
at present stage. However, the attempts concerning phonon 
properties in MNC lattices are less [11, 12, 14], compared 
with studies carried on momentum-conserving (MC) counter¬ 
parts exhibiting anomalous heat conduction [15-23]. More 
efforts should be devoted to the former. 

In a previous study [23], we have demonstrated the ad¬ 
vantages of a variational approach in describing renormalized 
phonons in MC nonlinear lattices with either symmetric or 
asymmetric nearest-neighbor potentials compared with exist¬ 


ing quasi-harmonic theories. Moreover, besides our method, 
we found that the nonlinear fluctuating hydrodynamics (NFH) 
[24] can also provide good estimation for the sound velocity 
of MC nonlinear lattices. Since the concept of sound veloc¬ 
ity is invalid in MNC nonlinear lattices, we deduce that the 
NFH can not be able to give information of the renormal¬ 
ized phonons, although it can still describe correlation func¬ 
tions of these systems [25]. Therefore, an interesting question 
arises whether the variational approach can be exploited for 
this class of models. 


The purpose of this work is thus twofold: (i)By putting the 
proposed variational approach onto a more general ground, 
we show that MNC nonlinear lattices can be tackled within 
the framework of the variational approach, which makes our 
theory a candidate for a unified phonon theory, and (ii) we 
determine temperature-dependent phonon properties of MNC 
nonlinear lattices using this method and confirm the validity of 
so-obtained results with the aids of molecule dynamics (MD) 
simulations as well as a tuning fork method [26]. As a demon¬ 
stration, we consider one of archetype ID nonlinear lattices in 
classical statistical mechanics, namely, ID o 4 lattice. We con¬ 
clude that the approach opens a new perspective to understand 
renormalized phonons in various nonlinear lattices. 


The paper is organized as follows. In Sec. II, we present 
the general aspects of the extended variational approach. Two 
intrinsic relation satisfied by the optimal reference systems 
are revealed such that general nonlinear lattice can be tack¬ 
led. Explicit forms of optimal variational parameters are also 
obtained. In Sec.Ill, we focus on a comprehensive study of a 
specific nonlinear lattice, namely, the ID lattice. We found 
a power-law temperature dependence of the phonon band gap 
valid in the weak stochasticity regime as well as a rigorous 
relation between ensemble averages of the nonlinear lattice. 
Finally, we briefly summarize our main findings in Sec.IV. 
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II. THE METHODOLOGIES and 


We consider ID nonlinear lattices described by the general 
Hamiltonian [27] 


N 


**= £ 


■77 + V(q n — <?n-i) + U(q n ) 


( 1 ) 


where N is the particle number, p n denotes the momentum of 
n-th particle, q n = x n — na denotes the displacement of n-th 
particle from its equilibrium position na with x n the abso¬ 
lute position and a the equilibrium distance for the interaction 
bond (for systems with zero pressure, it coincides with the lat¬ 
tice spacing), and V represents the interparticle potential, U 
is the on-site potential. We use periodic boundary conditions 
such that go = qN- For brevity and without loss of generality, 
we take m = 1 and a = 1 as the unit of mass and length, 
respectively. 

For simplicity, we introduce S n = q n — q n - i • The system 
evolves according to the equations of motion (EOMs) 


F > F 0 + (H - H 0 )p (7) 

as the variational principles, where F 0 is the Helmholtz free 
energy of the reference system, p 0 [= and p[= 

e — stanc j f or jjjg canonical measure of the reference 
system and the nonlinear lattice, respectively, /3t = 1 /T is 
the inverse temperature (we set ks = 1 ). 

The trial Hamiltonian Hq contains a set of adjustable pa¬ 
rameters {xi}, an optimal reference system can be obtained 
by varying those parameters such that bounds of the varia¬ 
tional principle go to a relative minimum or maximum. Gen¬ 
erally, without specific forms of H and Hq, the condition for 
the bounds [c.f., Eqs. ( 6 ) and (7)] to be stationary with respect 
to one of the parameters \i is equivalent to 

=°- (8) 

and 


Qn = Pn, (2) 

Pn = V'(6 n+1 )-V'(6 n )-U'(q n ), (3) 



(9) 


where the dot and the prime denote the time and spatial deriva¬ 
tive, respectively. Notice that lattices with asymmetric inter¬ 
particle potentials (U(<5) ^ V(—6)) will induce nonvanishing 
internal pressure, here we only consider lattices with symmet¬ 
ric ones. The extension to systems involving asymmetric po¬ 
tentials is straightforward [23]. Accordingly, nonlinear poten¬ 
tials with the following forms are of special interest: 


V(6„) 

OO p 

— /2s r2s 

~ ^ 2s n ’ 

S= 1 

(4) 

U (qn) 

OO 

_ V- 92k 2k 

- 2k qn ' 

(5) 


Coefficients are chosen in such a way that potentials exhibit a 
single-well character. 

In order to utilize a variational approach, the basic routine 
is to introduce a reference system with a trial Hamiltonian Hq 
and prepare the nonlinear system and the reference system 
at the same temperature, then determines optimal reference 
systems from variational principles [28], the resulting opti¬ 
mal variational parameters inevitably take temperature depen¬ 
dence. 

For systems with zero pressure, the Helmholtz free energy 
F determines its thermodynamic properties. Since nonlinear 
lattices in thermal equilibrium behave like weakly interacting 
renormalized phonons [18, 20 , 26], the free energy of the non¬ 
linear lattices can be divided into two parts, one is the contri¬ 
bution of the free renormalized phonon gas described by a har¬ 
monic Hamiltonian, the other comes from the phonon-phonon 
interaction. However, the interaction is complicated, approxi¬ 
mations should be taken. Intuitively, we adopt the well-known 
first-order cumulant inequalities of the free energy [29] 

( 6 ) 


for Eqs. ( 6 ) and (7), respectively, where (A; B) stands for 
the connected expectation value (cumulant): (AB) — ( A)(B}. 
These two identities result from intrinsic properties of the 
variational principles. Only those reference systems who sat¬ 
isfy one of the two identities can be regarded as optimal ref¬ 
erence systems. For later convenience, we refer the optimal 
system with the property Eq. ( 8 ) or Eq. (9) to the upper 
bound harmonic (UH) or lower bound harmonic (LH) system, 
respectively. 

The main concern of the present work is to obtain the 
phonon information in the nonlinear lattice from a theoreti¬ 
cal point of view. Phonon bears a solid basis only in harmonic 
systems, therefore, a harmonic reference system is chosen in 
the theory. Anharmonic references such as the Toda systems 
[30] are out of the present scope. For homogeneous nonlinear 
lattices with the chosen nearest-neighbor interactions as well 
as on-site potentials [c.f., Eqs. (4) and (5)], it is sufficient to 
consider the following trial harmonic Hamiltonian [31] 



in which fi 2 , and 7 are variational parameters with O 2 being 
the effective elastic constant, and 7 the strength of the on-site 
potential. Eqs. ( 8 ) and (9) determine their optimal values for 
the UH and LH system, respectively (see the following). The 
corresponding dispersion relation reads 

Wfe = 4D 2 sin 2 y+7 (11) 

with k is the wavenumber and a the lattice spacing. As can be 
seen, 7 also quantifies the phonon band gap. This dispersion 
relation with optimal parameters can be regarded as estima¬ 
tions of the actual dispersion in nonlinear lattices. 


F < Fq + (H — Hq) Po , 









3 


With the parameter set {xi} now being {H 2 , 7 }, using the 
identities Eqs. ( 8 ) and (9), we can obtain simple expressions 
for the corresponding optimal parameters (see the details in 
Appendix A). For the UH system, we will get the following 
self-consistent equations 


n 


2 

u 


7 u 


E hs(s 2 n s ) 

s=1 


Po 


(^)po 

E 92k (■<£*>, 

k= 1 


(&) t 


( 12 ) 

(13) 


With Q. 2 r and 7 u, estimations of the phonon dispersion [Eq. 
(11)] as well as phonon band gap via the UH system are thus 
determined. 

While for the LH system, we have similar forms as follows 


n 


2 

L 


7 L 


E f2s(6 2 n S ) 


P 



Eft* (<£*>„ 

k =1 

(qI) p 


(14) 

(15) 


However, unlike V.j- and 7 jj, these two parameters are to¬ 
tally determined by ensemble averages of the nonlinear lat¬ 
tice. With i 2 2 and 7 l, we can find another estimation of the 
phonon band gap and the phonon dispersion. It is worthwhile 
to mention that the above results can be applied to MC lat¬ 
tices, actually, fly and iij reduce to the known results for the 
FPU-/3 lattice [23], 


III. THE 4 > 4 LATTICE 

To verify the predictions of the variational approach, we 
consider the ID </> 4 lattice described by the following poten¬ 
tials [7-9] 


V(6) = ^5 2 , U(q) = jq\ (16) 

which corresponds to a special case of Eqs. (4) and (5) with 
/2 = K and <j,\ = A, respectively, and all other terms vanish. 
It is evident that only the on-site potential exhibits a nonlinear¬ 
ity. Being one of archetype ID nonlinear lattices in classical 
statistical mechanics, this MNC lattice has two advantages, 
on the one hand, the existence of renormalized phonons in 
this model has been confirmed by the use of the tuning fork 
method very recently [26], thus we can compare theoretical 
predictions with numerical results, on the other hand, the dy¬ 
namics of this model has been studied extensively, especially 
the existence of a strong stochasticity threshold (SST) which 
separates different phase-space structure as well as dynamic 
behaviors [32, 33]. It turns out that the temperature depen¬ 
dence of phonon band gap will manifest such transition be¬ 
havior. The SST also determines the validity regime of the 
optimal harmonic reference systems. In the following, we will 


give a detailed study of the 7 4 lattice by using our variational 
approach together with numerical simulations. For simplicity 
and without loss of generality, we only present details with 
K = 1 and A = 1, results with other parameter sets share 
similarities. 


A. Predictions of the variational approach 


Before going into details, we observe that f Ifj and O 2 [c.f., 
Eqs. (12) and (14)] should equal 1 by noting that the interpar¬ 
ticle potential of the ID </> 4 lattice is quadratic. Thus only 7 
needs further efforts. 

Firstly, we discuss the result obtained in the UH system. 
The numerator of Eq. (13) consists of only one term (c/ 4 ) pu . 
For a quadratic Hamiltonian, it is straightforward to show that 
(g 4 )po = 3(tj , 2 )p 0 [31], then the self-consistent equation [Eq. 
(13)] reduces to 

7 u = 3 (q 2 n ) P0 . (17) 

Note that (g 2 } Po can be directly computed from the Helmholtz 
free energy Fq via 

Z>n)po = 2 (18) 


For the harmonic system, F t) can be expressed as 
~Pt 1 E k l n EE with = 4 sin 2 + 7 , the above equa¬ 
tion thus leads to 


7 u 


3 x - 1 


(19) 


In the large N limit, the summation can be replaced by an 
integral, after some arrangements, the self-consistent equation 
[Eq. (19)] turns into a quartic equation 

7 ^+ 47 ^ = -^. ( 20 ) 

Pt 


The temperature dependence of 7 jj is explicit. 

Predictions from the LH system are quite straightforward. 
Eq. (15) implies that jl has the following form 


1L = 



( 21 ) 


which coincides with the prediction of the effective phonon 
theory (EPT) [19]. The temperature dependence of 7 ^ comes 
from temperature-dependent ensemble averages of the nonlin¬ 
ear lattice. 


B. Low temperature limit 

In the low temperature limit, by adopting the transfer in¬ 
tegral operator method (TIOM) [34—36], we find that(see the 
details in Appendix B) 

{q 2 ) p ~ 0.4477 T 2 / 3 , (22) 

(q 4 ) p ~ 0.5458 T 4/3 , (23) 
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which are consistent with results obtained from classical field 
theory [37]. We then deduce that 7 ^ has a following power- 
law temperature dependence 

7 L ~ 1.22T 2/3 . (24) 

The validity regime of this low temperature behavior will be 
identified by using MD simulations together with the SST in 
the following subsection. 

C. Numerical simulations 

In this part, we utilize MD simulations to verify the above 
predictions. A symplectic integrator SABA 2 with a correc¬ 
tor SABA 2 C [38] is adopted to integrate the EOMs of the </> 4 
lattice with a time step h = 0.02. We take an initial condi¬ 
tion such that the displacement of every particle is set to be 
zero and their velocities are randomly chosen from a Gaus¬ 
sian distribution at temperature T, after the initialization, a 
transient time of order 10 ' is used to equilibrate the system 
with N = 1024. 

The temperature dependence of 7 are shown in Fig. 1. We 
determine 7 jj by numerically solving Eq. (20) and getting 
its positive real solution. As for 7 ^, we insert MD results for 
(g 4 ) p and {q^) P obtained at temperature T into Eq. (21) to ob¬ 
tain its value at that temperature. The low temperature behav- 



T 


FIG. 1: (Color online)Temperature dependence of 7 for ID 0 4 lat¬ 
tice. The solid green line denotes the dependence of 7 u [Eq. (20)]. 
The dashed-dotted red line stands for the behavior of 7 l in the low 
temperature limit [Eq. (24)]. The blue circles are values of 7 l at 
temperature T using MD results for {q^)p and {qt) P at the same 
temperature. The black dashed line indicates the temperature corre¬ 
sponds to the strong stochasticity threshold. 

ior of 7 l [Eq. 24] is depicted in the figure as a dashed-dotted 
line. It’s clearly seen that such a low temperature behavior 
deviates from simulation results of 7 /_ when T > 0.6. This 
is because the temperature dependence of (q^) P and (q 4 ) p 
obtained from the TIOM [c.f., Eqs. (22) and (23)] are no 
longer valid when T > 0.6 compared with MD results for 
these two quantities as shown in Fig. 2. Noting that a temper¬ 
ature value 0.6 corresponds to a roughly average energy den¬ 


sity 0.5, which is just the SST found in the </> 4 lattice [32, 33]. 
Therefore the so-obtained low temperature behavior of 7 ^ is 
valid only in the weak stochasticity regime, once the system 
enters a strong chaotic regime, the simple power-law tempera¬ 
ture dependence of (q^) P and (g 4 ) p are not enough to capture 
the fast dynamics in the phase space. 



FIG. 2: (Color online)Ensemble averages for ID q !> 4 lattice. The red 
circles and maroon squares are MD results for {q^) p and (q 4 ) p , re¬ 
spectively. The solid green line and dashed blue line denote the corre¬ 
sponding low temperature behaviors. The black dashed line indicates 
the temperature corresponds to the strong stochasticity threshold. 

Using the predicted 7 u and 7 l [c.f., Eqs. (20) and (21)], 
we can check the validity of the dispersion Eq. (11) of opti¬ 
mal harmonic reference systems. We choose the tuning fork 
method introduced in [26] to obtain MD results for actual dis¬ 
persion of the ID 7 4 lattice. The comparisons between theo¬ 
retical predictions and MD results are shown in Fig. 3. The 
figure shows results for three different temperatures. One of 
the temperatures (T = 0.1) is well below the SST, imply¬ 
ing that harmonic references should be good approximations. 
From Fig. 3(a), we can see that this is indeed the case. The 
predictions of the LH system and the UH system both show 
a perfect agreement with MD results, since the values of 7 ^ 
and 7 u differ slightly at T = 0.1 as depicted in Fig. 1. As 
the temperature increases to 0.8 which is comparable to the 
SST, discrepancies between theoretical predictions and MD 
results appear near k = 0 and k = n as displayed in Fig. 3(b). 
But in the middle range of k, the two harmonic systems still 
can offer good estimations for the dispersion. The inset shows 
that the LH system is better than the UH system. When the 
temperature is well above the SST, the system is in a strong 
chaotic regime. Thus in Fig. 3(c), We find for T = 5 that 
our theoretical approximation Eq. ( 1 1) can not capture the ac¬ 
tual dispersion of renormalized phonon in the lattice. Noting 
the largest MFP in this case is almost 2 as shown in the inset, 
it means that the phonon-phonon interaction is very strong in 
this regime. Such that the first-order cumulant approximation 
adopting in the variational principle is no longer capable of 
dealing with lattices in the high temperature regime, higher- 
order corrections to the free energy should be taken into ac- 
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FIG. 3: (Color online)The dispersion for </> 4 lattice with length 
N = 1024 at different temperatures, (a): a comparison between 
MD results (red circles) and the approximation Eq. (11) either with 
7 l (solid green line) or •yu (dashed-dotted blue line) at T = 0.1. 
The inset shows details in the middle region of dispersion, (b): a 
comparison between MD results (red circles) and the approximation 
Eq. (11) either with 7 l (solid green line) or 7 u (dashed-dotted blue 
line) at T — 0.8. The inset shows details in the middle region of 
dispersion, (c): a comparison between MD results (red circles) and 
the approximation Eq. (11) either with 7 l (solid green line) or 7 u 
(dashed-dotted blue line) at T = 5. The inset shows MD results for 
mean-free-path (MFP) of the ID </> 4 lattice at T = 5. 


count. 


D. Exact relations for ensemble averages 

More surprisingly, although the variational approach takes 
an approximate way to get information of nonlinear lattices, 
we find that it can produce some rigorous results for the non¬ 


linear lattice. Noting the identity of the LH system Eq. (A14) 
and using Eq. (18), we find 77 also satisfies the following 
equation 

7! + 47i = ^I)|' ^ 

By combining Eq. (21), an identity of the ID (f> 4 lattice is 
revealed 

(q 4 n ) 2 p +Hq A n ) P (Q 2 n) P = T 2 . (26) 

Obviously, it is distinct from the simple relation in quadratic 
Hamiltonians, namely, (q 4 ) p = 3 {q 2 ) 2 - Moreover, it is worth¬ 
while to mention that this relation is valid in the whole tem¬ 
perature range even though (q^)p and {q 4 ) p have different 
temperature dependence in the weak and strong stochasticity 
regime separated by the SST. 

If we combine Eq. 26 with the generalized equipartition 
theorem [Eq. (A3)] for the </> 4 lattice with K = 1, A = 1 

(2 ql - 2q n q n - 1 ) p + (q 4 ) p = T, (27) 

we will further get relations hold between various ensemble 
averages, for instance, 

(qI) 2 p = ^(<7ngn-l)p + (q n qn-i)% (28) 

{Qn) pi^nQu—l) p = {Qn QnQn—l) p- (29) 

Therefore, in this model, it is apparent that the ensemble aver¬ 
ages of on-site quantities, e.g., q 2 , can be fully determined by 
the ensemble averages of non-on-site terms such as q n q n — i- 

IV. SUMMARY 

In summary, we have extended a previously proposed vari¬ 
ational approach to study phonon properties of general non¬ 
linear lattices with on-site potentials, which makes our theory 
a possible candidate for a unified phonon theory. Intrinsic re¬ 
lations characterizing optimal reference system, namely, the 
LH system and the UH system, are revealed, from which we 
can obtain explicit forms for optimal variational parameters. 
Estimations for the phonon bad gap as well as the phonon dis¬ 
persion are also obtained. 

As a specific case, we present a thorough study of the 
ID </> 4 lattice by utilizing the variational approach. In the 
low temperature regime, a power-law temperature dependence 
for 77 qualifying the phonon band gap is found by using a 
complementary method, namely, the transfer integral operator 
method. By compared with results obtained via MD simu¬ 
lations, we find that such low temperature behavior is valid 
only in the weak stochasticity regime where the temperature 
is below the SST. 

We further compare the theoretical predictions of the 
phonon dispersion with MD results for the actual phonon dis¬ 
persion of the ID <fi 4 lattice. At a low temperature below the 
SST, the agreements between theoretical and numerical re¬ 
sults are good. At an intermediate temperature which is com¬ 
parable to the SST, a discrepancy begins to appear between 
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theoretical and numerical results. We find that the LH system 
works better than the UH system. When the system enters 
a high temperature regime, or equivalently, a strong chaotic 
regime, our theoretical predictions fails, it seems that higher 
order corrections to the free energy must be taken into ac¬ 
count. 

More surprisingly, we find the variational approach can pro¬ 
duce rigorous results for the nonlinear lattice, such as an exact 
relation between ensemble averages of q^ and (/ 4 that is valid 
ranging from low temperature to high temperature regime re¬ 
gardless of the existence of SST. This finding together with 
the generalized equipartition theorem further lead to interest¬ 
ing relations hold between ensemble averages of the ID </> 4 
lattice. We then become aware of that ensemble averages of 
on-site quantities can be related to ones of non-on-site terms 
in the lattice. 
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Appendix A: Expressions of optimal parameters 

In this part, we derive explicit expressions of optimal vari¬ 
ational parameters for the UH and LH system using the iden¬ 
tities Eqs. ( 8 ) and (9), respectively. To do this, we should 
utilize the virial identity [39] 

(V • f) p = Pt (f • VH) p (Al) 

for any polynomial f which depends only on coordinates. If 
we take a special choice f = q with q = (91, 92 , ■'' > Qn) 
a vector consists of all the coordinates of the system with N 
particles, the above equality results in 

Pt (q- ViT) p = N (A2) 

or, specifically, 

= T ' (A3) 

which is the generalized equipartition theorem [40]. 

We firstly focus on the UH system. Choosing f = q{H — 
Ho) in Eq. (Al) with the canonical measure of the UH system, 
we get 

N(H-H 0 ) po + (q-V(H-H 0 )) po 
= Pt((H- H 0 )q • VHo) Po , (A4) 

On the other hand, we note that Eq. (A2) holds for the UH 
system 

(A5) 


Substituting this equation into Eq. (A4), we obtain 

<9 • V(fF - H 0 )) po = p T ((H- Ho ); q • VH 0 ) po • (A 6 ) 

The potential of the reference system is a polynomial function 
of the coordinates, then the scaling operation on Hq can be ex¬ 
pressed in terms of derivatives with respect to the variational 
parameters { \i }, i e., 

_ £jTT 

9-VJTo = X)Qi(x) 7 r 1 (A7) 

i ° Xi 

with a set of coefficients {Qi(x)} independent of coordinates, 
thus the right-hand-side of Eq. (A 6 ) vanishes due to Eq. ( 8 ) 
satisfied by the UH system, we are left with 

(<?• Vif) po = (q- V7T 0 ) PO ■ (A 8 ) 

With the concrete forms of H and Hq, the above identity gives 

Y hs ( S l S ) po + Y {q 2 n) po 

s—1 k— 1 

-^<^> P0 -7 u(q 2 n ) po = 0. (A9) 

Notice the facts that if H is a harmonic Hamiltonian, i.e., only 
/2 and (j ‘2 are nonvanishing in the potentials, the above equa¬ 
tion should immediately lead to O'/- and 7 u equal fi and < 72 , 
respectively, and also the nonlinearity of the on-site potential 
can not renormalize the interparticle coupling strength [31]. 
The above equation holds generally if and only if we let 

EMe> P0 =^<o P0 , (mo) 

S — 1 

E-92fe(9n fc ) po = lu{ql) po - (All) 

fc =1 

These two equations lead to expressions of 11/- and 7 u [c.f., 
Eqs. (12) and (13)], respectively. 

We then turn to the LH system. Applying the generalized 
equipartition theorem [Eq. (A3)] to the nonlinear system and 
the LH system, we obtain 

Y^( S n) P + Y92 k (q 2 n k ) p = T, (A12) 

. 9=1 k= 1 

nl{5l) po +lL(ql) po = T, (A13) 

respectively. For the LH system, the identity Eq. (9) with 
Xi = generates 

<£>* = <<a> (Ai4) 

<^} P0 = <«£>„■ ( A 15 ) 

Hence Eq. (A13) can be rewritten as 

tt|(O p + 7 L (q 2 n ) p =T. (A16) 

Subtracting Eq. (A16) from Eq. (A12), we get 

£Me> p -n 2 L (6l) p 

S= 1 

+ Y^(q 2k ) p -lL(q 2 n ) p = 0 . (A17) 

k= 1 


Pt (9 • VfT 0 ) Po = N. 
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Similarly, we should have 


where we have omitted the subscripts. 


E h.(sl’) r = a\(sl) p 

s—1 

(A18) 

(q 2 n) p = 1L (q 2 n) p ■ 

k=1 

(A19) 


These two equations produce Eqs. (14) and (15). 


Appendix B: Evaluation of ensemble averages with the TIOM 


In this part, we present in detail a study of the ID (f> i lattice 
via the TIOM. By adopting the TIOM, we should first intro¬ 
duce the following transfer integral equations [41] 

J dq n _ 1 e-P TK ^"-^'S’ l (q n -i) = 

J dqnMq^e-^n^’ 9 ^ =e~ 0Tei MQn-i), (B2) 


where K(q n , q n -i) satisfies 
^ ' K(.qni qn— l) 





T, and <b, are, respectively, the left-hand and right-hand 
normalized eigenfunctions corresponding to the eigenvalue 
e~d e i m if X takes a symmetric form, i.e., K(q n ,q n _ i) = 
K(q n _ i, q n ), we will have = \k*. So we choose 


E (q n , qn—i) 



(B4) 


In the continuum limit and low temperatures, the transfer 
integral Eq. (Bl) can be reduced to a pseudo-Schrodinger 
equation 


= E,'ii i (q), (B8) 

which describes a pure quartic oscillator with a temperature 
dependent coefficient. Introducing a transformation q = A q, 
the above equation becomes 




*<(Ag) 


X 2 Ei^i{Xq). (B9) 


We choose A 6 /3f, = 1 or A = pj. 1 ^ 3 such that the transformed 
quartic potential is temperature independent, which implies 

_i /q 

that the transformed ground state wave-function'T o(/3 t ' q) 
should also be temperature independent, thus we can denote 
it by tfio(q). Under the transformation, the averages [c.f. Eqs. 
(B6) and (B7)] become 


(<? 2 )p 

(<? 4 )p 


Pt 

Pt 


■2/3 J </>o(g)g 2 rfg 
/ 4>l(.q)dq 
4/3 f Poiq^q^dg 

J<t>o(q)dq 


Ci/3 T 2/3 , 

C2Pt A/ \ 


(BIO) 

(Bll) 


The orthogonality and completeness of 4', are guaranteed by 
the Sturm-Liouville theory for Fredholm integral equations 
with a symmetric kernel. Using this complete set, the ensem¬ 
ble average (g{q n )) P for a local quantity g(q n ) is equivalent 
to 

I *l(qn)g{qn)dq n 
/ (Qn)dq n 


in which ci and c-i are two temperature-independent coeffi¬ 
cients. 

Since the potential A q 4 is symmetric, the ground state 
wave-function should be an even function of q, we only need 
to know the positive q part of the function which approxi¬ 
mately reads [43] 


in the thermodynamic limit [31, 41, 42], where T'q is the 
eigenfunction corresponding to the ground state. So we get 


(<Z 2 )p 

(<7 4 )p 


J^(q)q 2 dq 
f*l(q)dq ■ 
S*Uq)q 4 dg 

m(Q)dq ’ 


p 0 (q) = exp 


(B6) 

(B7) where 

I 


y(x)dx 


(B12) 


. , V2 2 2 2 

y{x) = - —x + -e z 


.^ EIl/aW+I (3 )-fW 3, EVaW 


(B13) 


with z = ^x 3 and the exponential integral ElDz) defined by 


/ OO 

er zt t~^dt. 


(B14) 



















8 


Thus we can determine C\ and C 2 by numerically calculating 
the integrals with the above form of y(x), it turns out that the 


ci and C 2 equal 0.4477 and 0.5458, respectively. Therefore, 
we get the forms Eqs. (22) and (23) in our analysis. 
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